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By Wertheim method the exact solution of the Percus-Yevick integral 
equation for a system of particles with the "repulsive step potential" interacting 
( "collapsing" hard spheres) is obtained. On the basis of this solution the state 
equation for the "repulsive step potential" is built and determined, that the 
Percus-Yevick equation does not show the Van der Waalse loop for "collaps- 
ing" hard spheres. 

During the last 15 years the system of "collapsing" hard spheres (CHS) 
has been of the greatest interest of some experts in the field of physics of the 
condensed matter. The system of CHS is the system of particles with pair 
interaction potential 



where Vq is the positive constant, r is the particles range, R is the hard 
core diameter, a is the soft core diameter [see pQ and the references given 
there]. This interest is connected with a great application of the system of 
CHS in mathematical modeling isomorphic phase transitions, transforma- 
tions in colloid systems, in the study of curve of fusion, etc. [2]. By this time 
the calculations of thermodynamic characteristics of the system of CHS have 
been made by means of molecular dynamics and thermodynamic perturba- 
tion theory within the framework of spin-liquid approach [see for example [3] 
and the references given there]. 

However, there is another approach to the study of system with interaction 
in statistical mechanics, i.e. the solution of approximate integral equations 
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for the pair correlation function. We have not come across the case of so- 
lutions of integral equation for the pair distribution function of the system 
of CHS. In our opinion this research is necessary for the study of possibili- 
ties of integral equations method, because the use of this method in liquid 
theory for describing phase transitions is still the subject of many discus- 
sions [1]. In this work the system of CHS is being studied in the approach 
of Percus-Yevick (PY). It is an open secret that the PY equation admits 
the exact analytical solution in case of "hard spheres" potential [5]- [9]. Up 
to now it is the only analytical solution of nonlinear integral equation for 
the pair distribution function. On the basis of the method suggested in [6] 
we have obtained the exact analytical solution for a more complicated and 
"realistic" system of CHS for three-dimensional space when the number of 
particles and temperature are constant. Let's consider the PY equation in 
three diemensions (we follow the signs used in [6]) 

r(r) = 1 - p J r(r)/(r) dr + p J r(r')/ (r')r(r - r')e(r - r) dr' , (2) 

where e(x) = e~^ v ^ x \ f(x) = e(x) — 1, p is the density, and (3 = [kT)~ l . The 
pair distribution function g(x) and the direct correlation function C(x) of 
Ornstein and Zernike are in the PY approximation related to t(x) by 

g(x) = r(x)e(x), C(x) = r(x)f(x). (3) 

We make the equations dimensionless by defining x = (r/R), r\ = (nR 3 p/6), 
a' = a/R. Then we move to bipolar coordinates and integrate for angle vari- 
able for the repulsive step potential, and take the one side Laplace transform 
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of the equation satisfied by h(x) = xr(x). Defining 

l+a' 1 l+a' 

F(t) = - J xC(x)e~ tx dx = J h(x)e~ tx dx + e 2 J h(x)e~ tx dx, 
o o 1 

oo l+a' oo 

G(t) = J xg(x)e~ tx dx = e 1 J h(x)e~ tx dx + J h(x)e~ tx dx, 

1 1 l+a' 

l+a' 1 l+a' 

K = -F>(0) = - J x>C(x)dx = jx K x)dx + e 2 J xh(x)dx, 



(4) 



T(t) = J y(x)e~ tx dx, 
o 

l+a' l+a' 

y{x) = — y x'C(x')(x — x')g(x — x')dx' = ei€2 J h(x')h(x — x')dx' 

l+x 1+x 

ei = e-^°, e 2 = l-e a 
we obtain 

t [F(t) + G(t)} = 1 + 2 f 4vK + U V {[F(t) - F(-t)} G(t) - T(t) + T(-t)} . 

(5) 

Solving for G(t), this becomes 

1 + 24/7^ /N 12?7 r x % 

7^ F(i)-f[T(-t)-T(t)] 

G(t) = * ~ *- • (6) 

l + -l[F(-t)-F(t)] 

For further investigation we define the folowing function 

H(t) = t A G(t) { 1 + ~ F (-t) + ^ [T(-t) - T(t)]} + T(t) + T(-t). 

(7) 

Discussion like in [6] show that 

= Ax + A 2 t 2 + A 3 t 3 , 

where Ai, A 2 are constants. However, it is easy to check that H(t) is an even 
function, and so we define 

H(t) = A x + A 2 t 2 . 
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Not including G(t)F(—t) from (6) and (7) and turning the Laplace transform 
into the x < 1 + a', we gain the explicit expression for h(x) 



h(x) = C + C x x + C 2 x 2 + C 4 x 4 , 
where Cj is the solution of the following system of equations 

" ei /o (2) (/o (1) + e2 /o (2) ) - ±V* - - 2be ^o 2) 

&/ 3 (1) + be 2 fP - 126s/! + 6/ (1) + 3/ (2) (2e 2 + ae x 



C = b 
C\ — a 



C 2 = — b 
6 

Ca = — ab 
24 



and ei 



-pv 



, e 2 = 1 — ei,b = 12i]. The constants oe,fn\fn^ and y n as 



(8) 



(9) 



density, temperature and potential parametres Vq, R, R + a function can be 
obtained from 



f(i) 

J n 



x n h(x)dx 



l+o' 



= J x n h(x)dx 
1 

a = 1 + 2A V (/x^ + ea/^) 

a' 

Vn = J x n y(x)dx 



1+a' 



1+x 



(10) 



The system (9) must be completed with the condition C = 0, because 
h(0) = as it is followed from (2) and the nature of h(x). The system of 
"collapsing" hard spheres state equation can be shown as 



1 - 



pkT 6kT 



IT 



I™ 

(l + 4 V [ei/i(l) + e 2 (1 + a'f h(l + a')] ) 



dV 
dr 



(11) 



and inverse isothermic compressebility 

S) T w = 1 -»I c ^= x ^- < 12) 

It can be seen from (12) and (10) that if Vq > 

i.e. the Van der Waalse loop is absent in the isotherm. 
The gained solutions admits the generalization in case of n-step potential of 
particles. The calculations show that in case of n > 2 the Van der Waalse loop 
exists. Besides, the developed method allows to gain approximate solutions 
of the Percus-Yevick equation and, therefore, state equations for any smooth 
potential of particles. For this reason it is enough to appriximate the smooth 
potential with the n-step one. The systems of two-step potential are of the 
major interest nowadays [10]- [12]. The state equation of such systems in the 
Percus-Yevick approximation is the issue of our next article. 
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